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5_l [ Abstract 

We describe a combinatorial algorithm which, given a monotone and consistent sym- 
metric set function d on a finite set V in the sense of Rizzi RizOO , constructs a non 
trivial set S minimizing d(S, V \ S). This includes the possibility for the minimiza- 
tion of symmetric submodular functions. The presented algorithm requires at most as 
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much time as the one in |RizOO| . but depending on the function d, it may allow several 
improvements. 
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1 Introduction 



Let V be a finite set. A symmetric set function d is a map assigning a real number to each pair 
(yQ (S, T) of disjoint subsets of V satisfying d(S, T) = d(T, S). The function d is called monotone, 

if for every pair S, T of disjoint sets and T C T the inequality d(S*, T") < -i^, T) holds. It is 
called consistent, if for every triple R, S, T of pairwise disjoint set d(S, R) > d(T, R) implies 
d(S, BUT) > d(S U R, T). As the examples below indicate, the construction of a minimum 
bipartition of (V, d), ie. a nontrivial 1 subset S of V minimizing d(S, V \ S), is important for 
, many applications. 

ZS) Example 1.1. // G = (V, E) is a (hyper-)graph with weighted edges, the weight w(S,T) 

of edges between two disjoint sets of vertices is a monotone and consistent symmetric set 
function. A minimum bipartition of (V,d) is a minimum cut of the (hyper-)graph G. 'NI92a, 

Example 1.2. ///: 2 V — » R is submodular, ie. 

f(S) + f(T)>f(SUT) + f(SUT) 
for all subsets S,T C V , then the generalized connectivity function 

c f {S,T) :=f(S) + f(T)-f(SUT) 

is a monotone and consistent symmetric set function. If f is symmetric, ie. f(A) = f(V\A), 
then a minimum bipartition of Cf is a nontrivial set minimizing f. UJunSfA \Que9S\ QlizOOl 

A local version of the minimum bipartition problem, consists in the detection of a set S, 
such that d(S, V \ S) is minimal among all sets separating two given elements s and t, ie. 
s £ S and t S. We define 

\v,d) (s,t) = min {d(S, V\S)\seS,t^S} 



1 Nontrivial means ^ S ^ V 
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and 

\v,d) = nun{A(v j( Q(s,t) | s,t e V) . 

Hence, we want to find a set S, such that d(S, V\S) = \y,d) . 

Based on the work of Nagamochi and Ibaraki NI92b , NI92q] , Stoer and Wagner in |SW97| , 
and independently Frank in Fra94 , described an algorithm constructing a minimum cut of a 
weighted graph (exp. II. This algorithm was generalized by Queyranne |Que98| to the min- 
imization of symmetric submodular functions (exp. II. 2f) . and subsequently by Nagamochi and 
Ibaraki NI98 to a wider class of functions, satisfying a less restrictive symmetry condition. 

Another generalization was given by Rizzi in |R izOOj . He proved that the algorithm of 
Stoer/Wagner/Frank/Queyranne generalizes to the minimization of monotone and consistent 
symmetric set functions. Based on his work, we will describe a variation of the algorithm 
which provides possibilities to reduce the required time. Unfortunately, these improvements 
depend heavily on the function to be minimized. But we can guarantee, that the new 
algorithm requires at most the same time as that one described in RizOO . 

One example for the possible improvement can be found in |Bri05j . There, the author 
used the techniques presented in this paper to reduce the running time of the algorithm of 
Stoer and Wagner for the minimum cut of an integer weighted graph from 0(nm + n 2 logn) 
to 0(Sn 2 ), with 5 being the minimum degree 2 . 

We proceed as follows. First, we introduce lax-back orders, generalizing the max-back 
orders required by Rizzi in jRizOOj . Following that, we prove that the lax-back orders provide 
the possibility to identify several elements, hence reducing the size of V, without losing 
optimal sets. Then these results are used to describe and prove a new algorithm, bearing 
several possible improvements of the ones given in jRizOOj and |Que98| . 

2 Lax-Back Orders 

In the following, for each finite set V we write n = \V\ and a singleton set {v} C V will 
simply be denoted by v. 

Before we proceed, we make a simple and useful observation. 

Lemma 2.1. Let d be a symmetric set function on V and u, v, w elements in V and rel, 
then A(y )( j)(tt, v), \(y i( n(v,w) > r implies \(v,d)(u,w) > r. 

Proof. Let S be a subset of V with «gS and w g S, such that \fv,d) (u, w) — d(S, V \ S). 
If v e S, S separates v and w and hence d(S, V \ S) > X(vd)( v , w) > t. If v $ S we have 
d{S,V\S)>X {Vtd) (u,v)>r. □ 

Let V be a finite set and d a symmetric set function on V. An ordered pair (s, t) of 
elements of V is called good, if d(t, V \ t) — \i v<d \(s,t). As Rizzi proved in |RizOO| . such a 
pair (iVi_i, v n ) can be found by constructing a total order v\, . . . , v n on V, such that 

d(vi, . . . ,Vi-i}) > d(vj,{vi, . . .,Vi-i}) for 1 < i < j < n. 

An order of this type is called max-back order for (V,d). 

Lemma 2.2 (Rizzi [RizOO ). Let v±, . . . , v n be a max-back order for (V, d). Then (v n —i, v n ) 
is good for (V,d). 

Instead of using the original function d, we adapt it and introduce a threshold r. 

Lemma 2.3. Let d be a monotone and consistent symmetric set function and t a real number. 
Then d(S,T) = min {r, d(S, T)} is a monotone and consistent symmetric set function. 

2l Reduction' is a bit dangerous, since the running time is 0(Sn 2 ) and hence pseudo polynomial. Nonethe- 
less, if the edge weights are reasonable, one can expect the new algorithm to be faster than the old one. 



2 



Proof. Obviously d is symmetric. Now let S, T be two disjoint sets and T" C T. Then we 
have 

d(S, T') = min{r, d(S, T')} < min{r, d(S, T)} = d(S, T). 
Hence, d is monotone. 

Now let R, S, T be three pairwise disjoint sets such that d(S, R) > d(T, R). If d(S, R) < r, 
then d(T, R) < d(S, R). Due to the consistency of d, this leads to 

d(S, RUT) = min {r, d(S, R U T)} > min {r, d(S U i?, T)} = d(S U R,T). 

If on the other hand d(S, i?) > r, we have 

d(S, RUT) = min{r, d(5, i? U T)} > min{r, d(S, R)} = r. 

Furthermore, 

d(5 U i?, T) = min {r, d(S U R, T)} 

(t if d{SUR,T)>T 

~\d(SUR,T) if d(S U i?, T) < t 

< t < d(S,R\JT), 

concluding the proof of the consistency of d. □ 
A max-back order v\ , . . . , v n for (V, d) satisfies 

min{r, d(« i; . . . , > min{r, d(u,-, {wi, . . . for 1 < i < j < n. 

Such an order is called lax-back order with threshold t for (V,d). 

Corollary 2.4. Let d be a symmetric set function on V and v±, . . . , v n a lax-back order with 
threshold r for (V,d). Then the pair (v n -i,v n ) is r-good for (V,d), ie. for each set T with 
v n G T and v n -\ T , we have 

min{r, d({v n },V\ {v n })} = min {t, \( V ,d)(v n , v n -i)} ■ 

The proof of the following lemma is straightforward. 

Lemma 2.5. Let d be a monotone and consistent symmetric set function on V . 

1. Ifvi,..., v n is a max-back order for (V, d), then it is a lax-back order with threshold t 
for every r. 

2. Ifvi,... ,v n is a lax-back order for (V,d) with threshold t, then it is a lax-back order 
with threshold r' for any r' < r. 

As the preceeding lemma shows, we can interpret a max-back order as a lax back-order 
with threshold oo. 

3 Contraction and Lax-Back Orders 

In the following let d be a fixed monotone and consistent symmetric set function on V. In 
RizOO Rizzi used a max-back order for (V, d) to identify one good pair (v n -i,v n )- He then 
identified these two elements, obtaining a monotone and consistent symmetric set function 
d' on a smaller set V' . This process of identification or contraction may easily be extended 
to arbitrary partitions of V . 
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Let V = {Ui, . . . , Uk] be an arbitrary partition of V. For each element v £V : the class 
of v is denoted by [v], ie. [v] — Ui for v £ Ui. If, on the other hand, S is a set of classes of 
V, we write 

U5= (J [/.Cf 

for the set of members of classes in S. 

If d is a symmetric set function, we define the induced function d\> on V as 

d v (S,T) :=d(US, UT). 

It is easy to check, that (iy is monotone and consistent if d is. 

Lemma 3.1. Let d be a monotone and consistent symmetric set function on d. Ifvi, . . . , v n 
is a lax-back-order with threshold t for (V, d), then 

min{T, A(y;d)(uj-i,u,)} > min{r, d(v l} {v\, . . . ,^_i})} 

for 2 < i < n. 

Proof. Obviously, the restriction di of d to Vi :— {v\, . . .Vi} is a monotone and consistent 
symmetric set function on Vi and hence v± , . . . , Vi is a lax-back-order with threshold t for 
(Vi,di). By corollarv l2.4l we have 

min{r, d l (v l ,V i ^i)} = min {r, \rv u di)(vi-i, u,)} . 

Now let £> be an arbitrary subset of V - with Vi S S and ^ S". Then := SnVi separates 
Vi and Vi-i in and hence, due to the monotony of d, 

A (V r lA) (wi,«i_i) < di(5i, V$ \ 5i) = d(5i, V, \ .5^) < d(S, 1/ \ S), 

and thus 

A(y i ,d i )(^,«i-l) < \v,d)(vt,Vi_i). 
In combination with the observation made above, this leads to 

min{r, di(vi,Vi-i)} = min {r, A(v i ,d i )(^i-i, < {t, A(v,d)K, ■ 

□ 

Let Vi, . . . , v n be a lax-back order with threshold t for (V, d). If A(y^) < r, then no pair 
Vi) of elements with d(«i, Vi_i) > r can be separated by a set S with d(S, V\S) — \v,d)- 
Hence we may identify Vi and Vi—i, without increasing Xiv,d)- More precisely, we have the 
following result. 

Lemma 3.2. Let d be a monotone and consistent symmetric set function on V and r. IfV 
is a partition of V , satisfying Ary t< n(u,v) > r for each pair u,v of elements with [u] = [v], 
then for each pair s, t of elements of V with [s] ^ [t], we have 

min{r, A( V ,dv)( s >*)} = ™ n { T ,\v,d)(s,t)} , 

and 

min {r, X( V ,d v ) } = min {r, A (y , d) } . 

Proof. Let S be a subset of V with [s] S S and [t] ^ S. For 5 = US we have dv(£> V \ S) = 
d(S,V\S), and thus A( V ,d v )([s], [£]) > \v,d)(s,t), resulting in 

min {t, A(v,dv) } ^ min { r > A(v,<z) } • 
If \v,d)( s > t) > r, this implies 

min{r,A(v,<i v )([s], [i])} = min {r, A (y , d) (s, i)0} . 
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Algorithm 1: Lax-Back-Order 



Input: A consistent symmetric set function d on a finite set V, given by a lax oracle 

F and a threshold r 
Output: A lax-back order L = (v\, . . . , v n ) with threshold r for (V, d) 

Data: An ordered list L and a subset U of V 



1 L= («i); 

2 U = V\ui; 

3 while |[/ 1 > 1 do 



4 x := 0; 

5 forall it' G U do 

6 if V \ U; t) > t then 

7 L := (L, u), U :— U \ u; II Append u to L . 

8 x := r; 

9 else if \ U; r) > x then 
10 x:=F(u,V\U;t); 

n u := u'; 

12 end 

13 end 

14 if x < t then 

15 L := (L, u), U := U \u; II Append u to L . 

16 end 



17 end 



Now assume that \ry t( n(s, t) < t. In this case, no set S separating s and t with d(S, V \ 
S) = X(y,d)( s jt), can separate two elements u and v with [u] = [v], because \v,d) ( u ) v ) — T - 
Hence each partition [u] is contained in either S or V \ S. Therefore, S induces a subset S 
of V with [J<S = S and hence ^(v,d v )([s), [t]) < A(v,d)(s,i)j implying 

min{r,A( V ,d v )([s], [t])} = A( V ,d v )([s], [*]) <A(v,<j)(s,i) 

= min{r,A(y >rf )(s,£)} < min {r, A (v ^ v) ([s], [*])} 

The second part of the lemma is a direct consequence of the first part and the definition 
ofA(v,d)- □ 

Now let vi, . . . ,v n be a lax-back order with threshold r of (V,d). Then define V as the 
partition of V consisting of the classes of the transitive and symmetric closure of the relation 
Vi ~ Wj-i iff d(Vi, {vi, . . . ,Vi-i}) > r. 3 By lemma ETT1 V satisfies the condition of 13.21 and 
hence 

min {t, A( V ,d v ) } = min { T > \v,d) } • 
This simple observation is the main tool for the algorithm presented in the next section. 

4 Minimizing Symmetric Set Functions 

4.1 Construction of Lax-Back Orders 

Our algorithm for the calculation of a minimum bipartition of (V, d) , requires the construction 
of a lax-back order as a subroutine. To achieve this, we assume, that we have access to a 

' j Equivalently, we may say, that the classes of V are the maximal subsequences Vi,...,v^ such that 
d(vj, Vj) > r for i < j < k. 
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lax oracle F for d, ie. a program which returns min {r, d(S, T)} for any input of two disjoint 
subsets of V and a threshold r. 

The correctness of the algorithm is intuitively clear. Nevertheless we will prove it in 
detail. We show that, if an element u is added to L in lines 7 or 15, it satisfies the inequality 

min {r, d(u, V\U)}> min {r, V \ U)} 

for every u' 6 U. Since F \ U contains all elements already in L, this implies that L is an 
lax-back order with threshold r for (V, d) at the end of the algorithm. 
If u is appended to L in line 7, we have d(u, V \ U) > t and hence 

min {r, d(u, V \ U)} = r > min {r, d(u', F \ C7)} 

for each v! EU . 

If w is appended to L in line 15, no element was appended to L in line 7 in this round. 
Hence d(u, V \U) is obvious maximum among all u' G U . This concludes the proof of the 
correctnes of algorithm ^ 

In each round of the while loop, exactly |?7| calls of the oracle are required. Since in each 
round at least one element is removed from U, and since U begins with \V\ — 1 elements 
(line 2), at most "C" -1 ) calls to the oracle are made with n — \V\. Since it is possible, that 
more than one element is removed from U in a round, this bound may be very conservative, 
depending on the function d. As a result, a lax-back order can be calculated in time 0{n 2 Tp), 
where Tp is an upper bound for the execution time of the lax-oracle F. 

But, depending on the lax oracle, the runtime required for the construction of a lax-back 
order may be smaller. For example, as proven in jSW97j . the construction of a maximum 
adjacency order on a weighted, undirected graph G = {V,E), the analogue of a max-back 
order, requires time 0(m + nlogn), where n = \V\ and m = \E\. In |Bri05j this was 'reduced' 
to time 0(m+5n) for the first, and 0(8n) for each subsequent lax-back order with a threshold 
t < 5. Here S is the minimum degree of G. 

As the above example indicates, it is possible in many situations to calculate d(u, Vi) 
from d(u,Vi-\) much fater, than computing d{u, Vi) directly. In these cases a variation of 
algorithm H m the fashion of the algorithm described in Quc98 is more apropriate. 

In algorithm [3 we use a priority queue Q with threshold r, ie. Q contains value-key-pairs 
(v, k) 6 V x R and provides three operations. 

• The insert operation adds a pair (v, k) to Q. 

• The deljmax operation removes a pair (v, k) from Q such that 

k > min {r, max{fc' | (V, k') e Q}} . 

• The update-key operation updates the key k = d(u,Vi) to k = d{u 1 V^+i) if Vi was 
extracted from Q. 

The time required by algorithm [21 obviously depends on the type of thresholded priority 
queue used and the time required by update-key. But, since we are looking at the most general 
case, we simply assume that each lax-back order can be calculated in time Tlob G 0(n 2 Tp). 

4.2 Constructing a Minimum Bipartition 

Algorithm constructs a minimum bipartition, using algorithmnas a subroutine. Before we 
begin with the analysis of the algorithm, three remarks are appropriate: 

• In the given form, algorithm[3]first builds a max-back order, since r = oo. Alternatively, 
the initial threshold may be set to 

min {d{v, S\v) \ v G V } , 

resulting in a lax-back order, even in the first round. But in the worst case this requires 
n additional calls of the oracle. 
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Algorithm 2: Lax-Back-Order with Priority Queue 



Input: A consistent symmetric set function d on a finite set V, given by a lax oracle 

F and a threshold r 
Output: A lax-back order L = (v\, . . . , v n ) with threshold r for (V, d) 

Data: An ordered list L and a priority queue Q with threshold r. 

1 forall v E V do Q.insert(u, — oo); 

2 while \Q\ > do 

3 v := Q.deLmaxQ; 

4 L := [L,v)\ II Append u to L. 

5 forall u G Q do update_key (it); 

6 end 



Algorithm 3: OptimalSet 
Input: A consistent symmetric set function d on a finite set V 
Output: A subset S of V with d(S, V\S) = A (v ,d) 

Data: A partition V of ^ and a real number r 

1 S:=9: 

2 V := V"; 

3 t := oo; 

4 while |V| > 2 do 

5 ([vi], ■ ■ ■ , [vk]) '■= Lax-Back-Order(V, dv, r); 
e if d([v k ],V\[v k }) <rthen 

7 S:=[« fc ]; 

8 r := d([v k ], V\ [ufe]); 

9 end 

10 for i — 2, . . . ,k do 

n if d([ui], . . . , > t then 

12 Join [fj-i] and [uj]; 

13 end 

14 end 

15 end 



• Of course, we assume in lines 6, 8 and 11, that the d ([vi], {[i>i], • ■ ■ , [i^-i]}) were stored 
during the construction of the lax-back order. 

• We assume that lines 10-14 require O(n) time. 

Now we turn our attention to the correctness of algorithm |21 First, we prove that at the 
end of each execution of the body of the exterior loop (lines 5-14), we have 

\v.d) < r = d(S, V \ S) and Vu, V G V : [u] v = [v] v =>- X (V4) > r, (*) 

where [u]\> denotes the equivalence class of u G V in V. 

Observe that after the first execution of the body of the while loop, we have S = {v n } 
and t = d(v n , V \ v n ). Hence (*) is satisfied. 

Now assume that S, r and V are the 'values' before the execution of the loop and that 
they satisfy (*). If d([vk], V \ [vk]) > t, the values of S and r aren't changed and hence the 
first part of (*) is still valid at the end of the loop-body. If d([vk], V \ [v k ]) < t we obtain the 
new values r' = d([v k ],V\[v k ]) < t and 5' = [v k ]. Since d v ([v k ], V\ [v k ]) = d(S',V\S') = t', 
we obviously have X(v,d) ^ T> — d(S', V\S'). 
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In line 10-14, sequences [vi]v, . . . , [vj]v of classes in V are joined. This results in a partition 
V, such that each class [u]y of V is contained in a class of V". 

Now assume that [it]y = [v]y for two elements s,£ of V. If [u]\> = [i>]v> we have 
w) > r, since (*) was satisfied before the sets were joined. 

If [u]v 7^ [ v ]vj we assume that S* separates u and u with u E S and A(y,j)(w, u) = 
<f(5, V\S). If 5 also separates two elements s and £ with [s]v = [£]v, we have d(5, V\jSQ > r, 
by induction. 

If on the other hand, each class of V lies either in S or in V \ S, w.l.o.g. we have a 
sequence [u]v = [«i]v: • ■ ■ , [ w j]v = Mv of classes in V, such that 

d([vi]v, {[«i]v 5 • • • j [<>j-i]v}) > r for i < ^ < j. 
Therefore, lemma m*1 implies 

A(v,d v )(N]v, [wj_i]v) > t for i < Z < j. 

By lemma this implies 

\v,d)(vi,vi-x) > t for i < I < j 

and hence, by lemma TI~T\ \r Vtt n(u, v) > t, proving the validity of (*) at the end of the loop. 

Since r < dv,d v ([ w fc]vi V \ bfc]v) after line 9, at least the sets [«fc_i]v and [vk]v are joined 
in lines 10-14. Hence, the number of partitions in V decreases each round and after at most 
|V| — 1 rounds, we have |V| = 1, and the algorithm terminates. 

In this situation (*) implies X(y,d) < T — d(S, V \ S) and since [u] = [v] for all u, v G V, 
we have X(v,d)( u , v ) > T , leading to 

X (v ,d)=r = d(S,V\S). 

Theorem 4.1. Given a lax oracle F for a monotone and consistent symmetric set function 
d on V 

1. The algorithm^ returns a set SofV with d(S, V \ S) — \v,d) ^ n time 0(n 3 Tp), where 
Tp is an upper bound for the time required by the oracle F. 

2. If Tlbo is an upper bound for the construction of a lax-back order with threshold t < 
min {d(v, V \ v) | v G V} for (V,d), then a set S with d(S,V \ S) = \v,d)) can be 
constructed in time 0{n{Tp + Tlbo))- 

5 Conclusion 

We observed that for each monotone and consistent symmetric set function dona finite set 
V, the function d(S,T) = min {t, d(S, T)} is monotone and consistent, too. This fact was 
used to weaken the conditions on the oracle required for the construction of a set G S C V, 
minimizing d(S, V \ S). Instead of using a (strict) oracle F, we only required a lax oracle 
providing min{r, d(S, T)} for a given threshold r. This allowed several improvements: 

• Depending on d, a lax oracle may require less time than a strict oracle (an example can 
be found in |Bri05p . 

• The usage of lax-back orders instead of max-back orders and the fact that d is monotone, 
allows a possibly faster construction of the order, since the number of oracle calls per 
clement is reduced. 

• Since more than one identification is possible per round, the total number of rounds 
may be reduced. 

Unfortunately, all three improvements heavily depend on the function d. Hence, only a 
detailed analysis of special cases, may lead to a guaranteed improvement. In full generality, 
we can only guarantee, that the runtime of the algorithm presented in this paper is at least 
as fast as the one in RizOO]. 
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